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ABSTRACT 

Fatigue-crack growth data for 2024-T3 aluminum was used 
to evaluate the accuracy of Forman’s equation. Forman’s 
equation was found to predict from one-third to twice the 
observed cycles to a given crack length when applied to 
multiple cyclic stress ratios. An alternative approach 
was developed, based on the observed data, which is 
applicable to analog or digital fatigue-crack growth 
prediction. Data acquisition based on constant Increment 
of crack length has been found to yield Inadequate definition 
of the crack-growth curve In the low-cycle range for the 
successful application of the method. An expression was 
developed which extends the application of the crack-closure 
concept to the entire range of cyclic stress ratio. An 
effective root-mean-square stress was defined which shows 
Indications of good correlation with the magnitude of the 
crack-growth curves. 
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I. INTRODUCTION 



Modern aircraft, particularly military aircraft, must be 
designed for high performance and long service life. The 
attempt to achieve higher performance in terms of Increased 
payload, range, speed, or maneuverability implies a reduc- 
tion in gross aircraft weight, or in the structural fraction 
of that weight. This constraint results in structural 
components which carry higher loads per unit area and places 
greater demands on their structural integrity. 

The design values for the safe-life of non-lnspectable 
components and the inspection Intervals for other components 
are predicated upon the service life and the predicted load 
spectruiTi specified by the contracting agency. The aircraft 
components are sized to meet these design values using 
fatigue prediction techniques. The estimates of remaining 
life and maximum safe inspection Intervals of operational 
aircraft are frequently re-evaluated in light of the opera- 
tional load spectrum using similar methods. A factor of 
safety, or factor of uncertainty, is used in these calcula- 
tions by which the expected life of a given structural 
configuration is reduced. This factor of safety (PS) 
directly reflects the accuracy of the prediction method used. 
Large PS, indicating poor accuracy, results in excessive 
structural weight when compared to the optimum and penalises 
the aircraft performance throughout the service life. 
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optimistic estimates of the accuracy of a method lead 
ultimately to costly retrofits and reduced aircraft 
availability. The long term economic efficiency of an 
aircraft is therefore dependent upon the accuracy of the 
fatigue prediction methods available to the designer and 
operator; thus considerable motivation exists for the 
improvement of fatigue life prediction accuracy. 

The fall-safe design philosophy, which admits the 
existence of a crack in the structure, is one of the fatlgre 
prediction approaches in current use. Crack growth predic- 
tions, in units of cycles from an initial size to its lengj[?u 
when it precipitates catastrophic failure, are an essential) 
element in the fatigue analysis of today's airplanes. 
Inspection periods and rework intervals can be prescribed 
from this information. This thesis reports an investigatiia 
of data and theory used widely, if not exclusively, in 
making crack growth predictions for aluminum alloys. 
Specifically, Hudson's fatigue crack growth data [1] for 
202^-T3 aluminum alloy was examined in order to establish 
confidence limits on Forman's equation, which appears to bj- 
the most widely accepted crack growth prediction procedure 
used. Finally, an alternate method for expressing crack 
growth behavior is proposed to improve the accuracy of ci^ai: 
length predictions. 



i 
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II. BACKGROUND 



The fatigue life of aircraft components is difficult 
to predict due to the many factors Influencing the structural 
response to unsteady loading. These factors may be grouped 
into four categories; material characteristics, environmental 
effects, structural configuration, and service loading. 

The fatigue response to unsteady loads is statistical 
in nature and is dependent upon the time history of the 
loading and material response, where material response is 
defined to be cycles to a given crack length. The response 
history may be divided into three stages. The first stage, 
termed crack Initiation, is characterized by sub-microscopic 
cracks. This stage has been shown to be sensitive to a 
variety of microscopic damage mechanisms. The Intermediate 
stage, termed crack propagation, is characterized by the 
relatively slow growth of a macroscopic crack (a fev; grain 
diameters) under the Influence of the ordered external load. 
Pinal failure occupies a small fraction of the fatigue life 
of the material and a very high crack growth rate exists 
during this stage. 

A. MINER'S RULE OP CRACK INITIATION 

The obvious complexity of the problem has led investigators 
to make certain simplifying assumptions and/or approximations. 
Miner [2] proposed a linear damage accumulation over the crack 
Initiation stage of the fatigue life. This proposal defines 



13 




i 



f: . 

s 



I 








damage as the fraction of the Initiation life in cycles, N^, 

t h, 

expended at the i^ load level. An Increment of damage 

'h Vi 

would then take place with each cycle at the 1 load level 
according to l/N^^ or for n^ cycles, the damage would be 

Since the accumulation of damage was assumed to be 
linear. Miner’s Rule predicts the fatigue life of a 
component as 



L = Z n. 
1 ^ 



and damage as 



D = 



Z 

1 




( 1 ) 



( 2 ) 



At the completion of the crack Initiation stage. 



Z = 1 



( 3 ) 



Fatigue life data is presently available as vs 
curves where S Is the alternating component of stress and 
is the number of cycles to complete failure of the 
specimen rather than formation of a first crack. Application 
of Miner’s rule to random load tests have yielded damage 
predictions at specimen failure considerably different from 
unity. Some refinements to Miner's approach to damage 
prediction have taken the form of weighting functions on each 

mi 

term or variable powers on each term such as D = Z (nj^/Njj^) 



These modifications as well as competing damage hypotheses 
have ofttlmes resulted in Improved prediction accuracy for 
a given set of data; however, the Miner summation of damage 
continues to be the most widely accepted due to its simplicity 
and the minimal Improvement in accuracy offered by its more 
complicated alternatives. 

B . CRACK PROPAGATION 

After the crack is formed, fall-safe design practice of 
crack-growth predictions requires a knowledge of the crack 
growth rate and the load under which the cracked structure 
would fall, both of vjhlch are not available from any dam.age 
theory. Many investigators, as summarized in Ref. 3 , have 
attempted to correlate the crack growth per cycle to the 
gross section stress and the crack length, leading to 
propagation models of the form, 

da/dN = f(S,a,X) (4) 

where 2a is the total crack length, S is the gross section 
stress and X is a correction factor which compensates for 
finite width, rivet forces, etc. Paris and Erdogan [4] 
concluded that the stress Intensity factor (SIP) , which is 
a function of the same variables, was the parameter which 
best correlated with the crack growth rate observed. 
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C. STRESS INTENSITY FACTOR 

Mathematically, the stress Intensity factors are undeter- 
mined coefficients arising from the solution of the field 
equations for any cracked geometry. The normal stress at 
the crack tip, under moderate static loading, is given by 
this solution as 






^2 

g 

r 



ij 



(e) + 






h 



h,j(6) 



( 5 ) 



where K^, K 2 , and are the opening mode, in-plane shear 
mode, and out-of-plane shear mode SIP respectively, r is 
the radial coordinate referenced to the crack tip and 6 is 
the angular coordinate measured counter-clockwise from the 
advancing crack tip. The stress intensity factors may 
therefore be interpreted as parameters that reflect the 
redistribution of the stress field in a body resulting from 
the introduction of a crack. The magnitude of the stress 
Intensity factor is dependent upon the geometry of the body 
containing the crack, the size and location of the crack, 
and distribution and magnitude of the external loads on the 
body. The configuration most readily treated is the sheet 
specimen containing a central crack of length 2a under 
Mode 1 loading, for which 
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where is a width correction factor. The criterion for 
failure in the presence of a crack-like defect is that crack 
growth to failure (instability) will occur whenever the 
applied stresses, S, as given by equation (6) exceed some 
critical condition specified by K, = K . Paris and Ergodan 

-L C X 

concluded from the preceding argument that this condition 
may be extended to fatigue and applied to the slow growth 
phase. 

D. FORMAN’S EQUATION 

Forman [5] proposed a model for crack propagation which 
utilized the range of the SIF as suggested by Paris [6], 



AK 



K — K 
max ~ min 



(7) 



to account for the cyclic loading typical of the fatigue 
environment. This model included the instability condition 
for failure as a singularity such that the growth rate 
approaches Infinity as K approaches K . Forman’s 

IT13. X cr 

equation is written as 



da _ C(AK)^ 
dN (l-R)K^^ - AK 

where C and n were assumed to be undetermined material 
constants and R is the ratio of minimum stress to maximum 
stress (loading ratio). 
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E. HUDSON'S DATA 



Hudson [1] conducted a most extensive study of the 
material response of two widely used aircraft alloys, 

2024-T3 and 7075-T6 aluminum, to evaluate Forman's equation 
when applied to a large number of maximum stress conditions 
and loading ratio combinations. (The data obtained in 
Hudson's investigation was used as the data base for this 
study. ) 

The specimens used in Hudson's study were nominally 
12 inches in width, 35 Inches in length, and 0.090 inches in 
thickness. A notch 0.10 inches in length by 0.01 Inches in 
width was cut into the center of each specimen by an elec- 
trical discharge machining technique. A reference grid of 
0.050 inch spacing was photographically printed on the 
surface of each specimen to mark intervals in the path of 
the crack. 

"Fatigue-crack growth was observed through 10-power 
microscopes while Illuminating the specimen with stroboscopic 
light. The number of cycles required to propagate the crack 
to each grid line was recorded so that crack-propagation 
rates could be determined. Approximately two-thirds of the 
crack-propagation tests were conducted to failure. The 
remaining one-third were stopped before failure, and the 
cracked specimens were used in residual-static strength tests. 



National Aeronautics and Space Administration Technical 
Note D-5390, Effect of Stress Ratio on Fatlgue-Ci~‘ack Growth 
in 7075-T6 and 2024-T3 Alumin um Alloy Specimens , by C. 
Michael Hudson, p. 5, August 1969- 
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The results of Hudson's fe.tigue-crack-growth tests on 
the 2024-T3 specimens are presented In Table I (reproduced 
from Ref. 1). Hudson obtained the fatigue-crack propagation 
rates graphically from the crack-growth curves defined In 
Table I. The crack-growth curves reproduced from Ref. 1 
and shown as Figure 1 are typical. The results were 
presented as plots of the common logarithms of the experi- 
mental value of the crack-propagation rate as functions of 
the range of the SIF. Three models were then fitted to the 
data by least-squares techniques and compared to the data. 
Figure 2, from Ref. 1, is typical of the method used for the 
presentation of data and the evaluation of the quality of 
fit of empirical equations. Hudson concluded that Forman's 
equation produced an excellent fit for both materials whereas 
the other equations evaluated "...produced good correlation 
with the test data except at the high growth rates for the 
7075-T6 alloy" on the basis of graphical representations of 
which Figure 2 is typical. This conclusion at first appears 
to be justified on the basis of plots such as Figure 2, 
hovjever, the fatigue analyst is concerned with the accuracy 
with which a crack-propagation model predicts the number of 
cycles required for the crack to attain a given length. This 
is the basis upon which the fit should be evaluated. 
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Ill . DEVELOPMENT 



A. EVALUATION OF FORMAN’S EQUATION 

Hudson used the tangent form of the width correction 
factor. 



in evaluating the material constants, C and n, in Forman's 
equation. The differential equation resulting from the use 
of this form of width correction in Forman's equation may 
be Integrated (Appendix A) using the small angle approximation, 
tan 0 9 . This approximation is accurate to within three 

percent for crack lengths of less than 0.8 inches for 
Hudson's specimens. 

The predicted Interval of cycles between crack lengths, 

ANp, was calculated using Hudson's fit of Forman's equation 
and the Integrated form of Forman's equation becomes 




( 9 ) 








( 10 ) 



where 



(1-R) 



C 



1 



C[(1-R)S 



] 



n 



max 
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and AN is the predicted number of cycles required to 



propagate at half-crack of length a^ to length a. The ratio 
of the predicted to the observed cycle Increment is presented 
in Table II, where gross discrepancies exist for almost every 
case. 

The load cases were also fitted individually in order to 
determine the cause of the poor predictions reported in 
Table II. The expression for AK used in this analysis was 



where Brown and Srawley's [ 7 ] finite width correction for a 
centrally cracked sheet. 



was chosen on the basis of its reported accuracy. The test 
data for each load case were represented by a power series 
of the form 



which yields a very good approximation of the observed data 
as indicated by Figure 3* An analytic crack-growth rate 
expression v;as obtained by performing a differentiation of 
equation ( 13 ) with respect to cycles leading to 




( 11 ) 



X„ - 1.77 + .277(^) - .51(^)^ + 2.7(^)^ (12) 



a = p,N + p„N + . . . + p N ^ 
1 2 m 



(13) 
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+ 



(14) 



dN 



e 

p^e^N 




ej-l 



p e N 
m m 



e -1 
m 



Equations (11) and (13) v;ere applied to each load case and 
a least-squares fit of Forman’s equation to the data was 
performed. Figure 4 Is typical of the Individual fits thus 
obtained. 

A regular, periodic variation of the calculated crack- 
growth rate about the fitted line was observed which Is 
believed to Indicate that Forman's equation m.ay be an 
Incorrect or Incomplete functional form. These periodic 
variations appear as scatter when Forman's equation Is called 
to fit multiple load cases. An apparent lower bound of 
crack-growth rate may also be observed In Figure 4. The 
failure of Forman's equation to represent this threshold 
level probably accounts for the largest portion of prediction 
error. Just as for the Integration performed for Hudson's 
fit, consistent prediction errors were noted for the fits 
of the Individual load cases but of smaller magnitude. The 
periodic variation of the data about the fitted line and the 
apparent lower threshold of crack-growth rate led to the 
development of a new material response model. Equation (13) 
was not considered to be a pralctlcal alternative to Forman's 
equation In view of the explicit appearance of the variable 
N, which Is unknown In the practical application of any 
prediction model. Series reversion was likewise discounted 
due to non-integer exponents resulting from the fit. 
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B. MATERIAL RESPONSE FUNCTION 

A material response function was developed in which the 
response was assumed to be a piecewise continuous function 
of time, loading conditions, and crack geometry. The 
function was also assumed to possess piecewise continuous 
time derivatives. This assumption was based on the satis- 
factory representation of the data achieved by the power 
series approximation of the fatigue data, equation (13) j 
which has an infinite number of continuous time derivatives. 

Previous investigators have assumed that the fatigue- 
crack propagation rate, da/dN, has the functional form 

da/dN = f(S(N),a(N),X(N)) (15) 

where the variables have the meanings previously defined and 
are functions of time or, alternatively, cycles. Although 
the relationship between the independent variables N and t 
is linear, as expressed by 

N = wt ( 16 ) 

where w is the circular frequency of loading, a distinction 
in the growth rates da/dt and da/dN exists. The crack-growth 
rate da/dt is an instantaneous rate v/hich must be integrated 
over the period of a cycle in order to obtain the time 
average rate, da/dN. 
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Paris C6] concluded that the crack growth rate was best 
represented as a function of the range SIP as given by 
equation (11). The desired function must predict crack 
growth rate under constant amplitude and random loading with 
equal precision since constant amplitude loading is a special 
case of random loading. The stress range is therefore a 
constant for either load case taken on a cycle-by-cycle basis. 
This argument applies equally well to any other choice of 
the characteristic stress to be used in calculating the SIP. 
The expression for crack-growth rate could then be represented 
equally well by the functional form 

(17) 

( 18 ) 



since X is a function of crack geometry only. 

Introduction of the requirement for higher order deri- 
vatives can be satisfied by a linear summation of the assumed 
form 



Or 

§ = f2(^(N)) 

where 

X = X^(a,W)'\/a(N) 



= c 



dx , 

jTT + c< 



1 dN 



dfx 

dN^ 



+ 
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X 



• • • 



(19) 



where the = c^(S). Equation (19) Is a linear differential 
equation In N with constant coefficients and thus has a 
solution of exponential functions. Paris [6] noted that the 
crack growth rate Is a monotonlcally Increasing function 
of the stress Intensity factor. The crack length, a. 
Increases monotonlcally with time and the width correction, 

X, Increases monotonlcally with crack length, therefore the 
rate of change of the variable", x, with respect to cycles 
(time) must Increase monotonlcally with x. This type of 
response Is typical of positive feedback systems, which 
Implies that the roots of the characteristic equation 
obtained from equation (19) are real. A mathematical model 
of the desired material response function Is therefore the 
Prony series 



which satisfies all requirements heretofore set forth. 

This form of the material response function may be recast 
to the differential 



when the functional form of the c^^ Is determined. This 
differential equation would find utility If programmed for 
analog or hybrid computation using non-llnear elements 



X 



2 ^21-1 6xp[b2j^N] 



( 20 ) 





( 21 ) 
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eliminating the excessive digital computer time requirements 
typical of random loading fatigue calculations; however, it 
Is beyond the scope of this thesis. This approach to fatigue 
prediction calculations would also allow the use of discretized 
flight load data without Intermediate processing. 

The data sampling frequency used by Hudson, based on 
equal Increments of crack length, provided insufficient data 
in the low cycle range of the specimen lives. The crack 
growth rate progresses from a nearly constant value at small 
crack lengths to much larger values within a small Increment 
of crack length requiring a more frequent data sampling 
frequency in order to achieve a successful analysis of the 
exponential form of the fatigue response function through 
direct analysis. 

C. EFFECTIVE STRESS 

The extension of equation (20) to account for random load 
sequences requires that the constants be a function of the 
loading conditions. The representation of the piecewise 
continuous response function, x'(t), as the cycle average 
function x(N) leads to the representation of the piecewise 
continuous crack-growth rate, da(t)/dt, as the time average 
crack-growth rate, da(N)/dN. Consistency requires that the 
loading also be represented as a time average variable. 

Paris [l6] suggested the utility of the Root-Mean-Square 
(rms) gross section stress, as defined as 
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s 2 i 

rms T 



( 22 ) 



T 

/ S^(t) dt 

0 

where S(t) = Sm + Sa sin wt 

and Sm is the mean gross section stress during the cycle, 

Sa is the amplitude of the gross section varying stress 
during the cycle, w is the circular frequency and T is the 
period of the cycle. 

The entire magnitude of the stress described above does 
not contribute to crack propagation. Elber [8] has noted 
that plastic yielding at the crack tip results in a compres- 
sive residual stress upon unloading. This compressive stress 
at the crack tip delays the opening of the crack. Elber 
defined the crack closure parameter, U, as 

u = ~ 3°P (23) 

^max min 

where S^p is the gross section stress at which the crack 
opens. He found that U seemed to be a linear function of 
loading ratio for 2024-T3 aluminum as represented by 

U = 0.4R + 0.5 (2^1) 

Since no crack propagation occurs while the crack is 
closed, an effective rms gross section stress may be defined 
as 
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where 



rms 



1 

/ 

0 



2 ^ 
e T 



+ S.sin 
m a 



2-nt 

T 




( 25 ) 



subject to Sg ^ 0 . 



The nature of the sinusoidal function is such that the limits 
of integration can be shifted to write (25) more conveniently 
in terms of the damage effective portion of the cycle as 
illustrated in Figure 5. This leads to 

1 ,J_ 1 _L 

T 2T T"2T 



rms. 



=/ 



c, 2 dt 

^e T 



/ 



0 ^ 
T 



( 26 ) 



1 _ _f_ 
¥ 2t 



1 +X 
ir^2t 



where f/T is the fraction of the cycle during which the 
crack is open and is defined by (See Appendix B) 



f/T = 



cos ^(1-2U) 
7T 



(27) 



Elber’s equation for U(R) leads to effective cycle 
fractions v;hich are undefined for load ratios less than -1.25. 
This unnatural behavior in the neighborhood of a subset of 
Hudson’s data motivated a search for a more satisfactory 
representation of the crack closure phenomenon. Figure 7 
is a schematic of the loading conditions for increasingly 
negative load ratios. Reasoning from this figure indicates 
that, in the limit, as S approaches zero (R approaches 
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negative infinity) , the fraction of the period over which 
the crack is open approaches zero. This physical constraint 
is satisfied by an equation of the form 



f ^ A 
T B - R 



( 28 ) 



where A and B were assumed to be undetermined material 
constants. Substitution of equation (28) into the defini- 
tion of f/T, equation (27), leads to 



U = I [1 - cos TT 



(29) 



Elber’s data [8] was fitted to the function for the 
effective cycle fraction, equation (28), as shown in Figure 
6. The form satisfies the physical constraint . and yields 
a smaller x error than equation (24) when compared in 
Figure 7. Figure 7 is a comparison of the functional 
representation of U(R) as developed in this study with that 
of Elber. Both forms appear to be incomplete in that neither 
form considers the effects of absolute load level or crack 
geometry on the crack closure phenomenon. 
calculated on the basis of equations (25) and (28). 

D. HYPERBOLIC REPRESENTATION OP FATIGUE -CRACK GROWTH DATA 
On the basis of the arguments presented with regard to 
the exponential form of the material response function and 
the observed similarity of the exponential and hyperbolic 
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curve forms, Hudson’s fatigue-crack growth data was fitted 



to the general form of the second order equation. 



Ax^ + BxN + CN^ + Dx + EN 



1 



(30) 



Solving for N, Table III is a listing of the coefficients of 
the equation 



ordered with decreasing , vjhich resulted from the fit 

of equation ( 30 ). 

Figure 8 indicates the quality typical of these coeffi- 
cients. The variation in sign which occurs in Table III 
reflects both the translation of the x axis (or the N 
ordinate of the center of the fitted hyperbola) and the 
frequent occurence of parabolic, rather than hyperbolic, 
coefficients. The parabolic coefficients (sign of A = sign 
of C in equation (30)) resulted from the paucity of data in 
the low-cycle range of the specimen fatigue life. The 
appearance of parabolic coefficients is unacceptable from 
a physical standpoint as the parabolic form implies that 
cycles to a given crack length is not a single valued 
function, as shown schematically in Figure 9a. The trans- 
lation of the center of the assumed hyperbola was an 
unavoidable result of the presence of the original machined 
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(31) 
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notch. Figure 9b is a schematic representation of the 
effect of load level and the presence of the artificial 
notch on the translation of the x axis. 



E. NORMALIZATION OP THE DATA 

A unified fatigue prediction method may be formulated 

by combining the capabilities of the damage accumulation 

method in estimating crack initiation with the fracture 

mechanics approach applied over the crack propagation and 

failure stages. This usage demands the normalization of 

all load cases to a common reference crack length which 

should be the minimum detectable crack length defining 

initiation. This normalization vms accomplished for each 

(hyperbolic) load case by estimating the cyclic interval 

betv/een Hudson's reference point, typically 0.10 inch, and 

the assumed Initiation points, 0.003 inch. The estimate was 

made on the basis of equation (31) using the d^^ previously 

calculated and was then applied as an increment to the 

observed cycles to a given crack length for an unnotched 

specimen. The modified data sets were again fit to equation 

(30) yielding estimated coefficients, d^*, of the unnotched 

sheet which are hypothesized to correlate with „ . 

rmse 

Figures 10 through 1^ indicate the degree of this correlation 
and suggest that the appropriate functional relation may be 



di* = 



ttj S 
1 rms. 



(32) 
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An analy.tlc expression for the crack-growth rate, da/dN, 



may be obtained by the chain differentiation 



da _ ^ ^ _ 1 1 

dN dx dN ('^'1 /dN^ 
Ma'' Mx'' 



(33) 



where 



dx _ 1 / ^ 4 . „ 

31 • — ‘ 2 ^ ® di’ 

a 



( 34 ) 



and 



da 



4 (2doX + di.) 



d. - 






d^x + dj^x + d^ 



(35) 



This expression for da/dN approaches the constant (threshold) 
crack-growth rate observed in Figure 4 at small crack lengths 
and satisfies the instability condition given that equation 
(34) has a zero. 

The zero oT equation (34) yields x^^^^ and leads to 



K,. o 
ins 



= S „ X. 
max ins 



(36) 



where the instability condition has been assumed to take 

place at the instant in the last cycle at which the load 

has attained its maximum value. Figure 15 depicts the 

variation of the ratio K.^^/K . This apparent reduction of 

Ins cr 

the level of the stress intensity factor for unstable crack 
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propagation Is attributed to dynamic effects on propagation 
and has the effect of minimizing the Influence of the 
singularity condition in Forman’s equation. 

The data fits to equation (30) yielded no real solution 
to equation (3^) for some load cases. These load cases 
appeared to be the dynamic tests which were terminated well 
before failure and therefore contained insufficient 
information regarding the shape of the response function in 
the hlgh-cycle region to give the correct physical response. 
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IV. CONCLUSIONS 



Hudson’s [1] fatigue-crack growth data for 2024-T3 
aluminum was used to evaluate the accuracy of Forman's 
equation. An alternative approach to the representation 
of the propagation process was developed and analyzed. An 
expression was developed which extends Fiber's crack-closure 
concept to the possible load ratio range. The following 
conclusions can be drawn from this study. 

1. Forman's equation does not appear to satisfactorily 
predict crack propagation in 2024-T3 when evaluated on the 
basis of the observed data. 

2. Fiber's crack-closure concept was successfully 
extended to the entire load ratio range in 2024-T3 by the 
relation 



U 



sin 



/ .526 

'^2.169 - R 



) 



3. The development of a crack-growth rate expression 
based on empirical analysis of the observed data met with 
limited success due to the coarse data resolution in the 
low-cycle range. 

4. The stress intensity factor at crack instability 
under dynamic loading conditions was observed to decrease 
with decreasing values of effective rms gross-section 
stress, asymptotically approaching a value of approximately 
one-third the value measured under static loading. 
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IV. RECOMMENDATIONS FOR FUTURE RESEARCH 



The material response of 2024-T3 aluminum has been shown 
to be capable of representation by a continuous response 
function. The accumulation of a number of test sets with 
more frequent data sampling In the low-cycle phase of the 
crack propagation life would allow a determination of a 
material response model for which statistical reliability 
limits could be determined. Any such model must be based 
upon the observed data, as compared to derived data, for 
the reliability limits to be meaningful. 

The development of the exponential response model will 
allow the determination of the governing differential equation 
for crack propagarion under dynamic loading by the application 
of parameter Identification techniques. Such a development 
would allow the calculation of the expected crack length 
under random loading by the use of analog, rather than 
digital, computation devices. 
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APPENDIX A. INTEGRATION OF FORMAN'S EQUATION 



Forman's equation is written 



^ = C(AK)^ 



dN (l-R)K^j, - AK 



(A-1) 



where the stress intensity factor, K, is defined as 



and the range of the stress intensity factor is written 



may be used, allowing the small angle approximation. This 
approximation, tan 6=0, yields errors of less than three 
percent for ratios of crack half-length to width of less 
than 1/15. 

Limiting the integration to this range allows AK to be 
expressed as 



K ^ SX-\f^ 



(A-2) 



AK = ASX^,^^ra' 



(A-3) 



The tangent form of the width correction factor 




(A-il) 



AK = AS^ya* 



(A-5) 
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The cyclic stress ratio, R, is defined 



R = 



max 

^mln 



(A-6) 



Therefore AS = (l-R)S 



max 



(A-7) 



Forman’s equation can be rewritten as 



^ ^max ^ 

(l-R)K^^ - S a^^^) 



(A-8) 



cr max 



for the limitations Imposed. 

Treating this expression as a differential as proposed 
by Forman [5] 



dN = 



(1-R)(K - S a^^^) 

cr max 



C(l-R)^S „ ^a^^^ 



da 



(A-9) 



max 



which may be integrated between crack half-lengths a^ and a 



a K - S a^/2 

r da 



P C(1-R)^S 



max 



W2 



(A-10) 



Carrying out the indicated integration yields 



AN„ = 



(1-R) 






2-n 



a 

cr 



2-n 



3-n _a 



S . a 
max 

3-n 



(A-11) 
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Simplification of equation (A-11) leads to 





2-n 




K S -V 
cr maxV o 

n-2 n-3 



2~n „ 

o 2 ^cr 



max 






n-3 



(A-12) 



where 



^1 = 



(1-R) 
C(1-R)^S 



n 



max 



(A-13) 



and an Is the cyclic increment predicted to elapse during 
the extension of a crack of half-length a^ to a half-length 
of a when subjected to constant amplitude loading conditions 
R and S 

max 
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APPENDIX B. THE DERIVATION OF THE RELATIONSHIP BETWEEN 
THE DAMAGE EFFECTIVE CYCLE FRACTION AND 
THE CRACK-CLOSURE PARAMETER 

Crack propagation occurs only when the fatigue crack 
is open, allowing the definition of a damage effective 
portion of a cycle (period) . This fraction of the period 
is bounded by the limits t^/T and t 2 /T representing the 
instants of crack opening and crack closing, respectively. 
The damage effective fraction is then 

f/T = t2/T - t^/T (B-1) 

These limits arc symmetric about the maximum load instant 
for constant amplitude loading. Assuming a loading function 
of the form 

S(t/T) = S + S^ sin (2TTt/T) (B-2) 

m cL 

the limits of the damage effective fraction of the cycle 
may be written 

t^/T = 1/H - f/(2T) (B-3a) 



and 



t2/T = 1/4 + f/(2T) 



(B-3b) 



The crack opens during the loading portion of the period 
when the effective gross section stress is zero. The crack- 
closure concept modifies this stress by the inclusion of 
the residual stresses resulting from local tensile yielding 
at the crack tip. Defining S^p to be the gross section 
stress at which the crack opens leads to 

Sg (t/T) = - S^p + sin (27rt/T) (B-4) 

The Instant of crack opening is therefore defined by 

Se (t^/T) = - S^p + sin (2Trt^/T) = 0 (B-5) 



or 



s — s 

(-^2_ 21) = sin [2 tt( 1/4 - f/2T)] 

■^a 



(B-6) 



Trigonometric reduction yields 



S - S 

(_OP^ 21) = cos (TTf/T) 



(B-7) 



and the damage effective fraction of the cycle is defined 
in terms of the loading and closure parameters 



f/T ^ cos ^/tt 



i|0 



(B-8) 



Rewriting Equation (B-8) in terms of S „ and S . 

max min 



yields 



f/T = 



cos ^/ir 



S + S 

<j / m ax mi n \ 

^op ~ ^ 2 ^ 

(^max ~ ^mln) 



(B-9) 



which simplifies to 



f/T = COS~^/lT 



1 - 



2(S - S ) 

max op 

S — S 
max min 



(B-10) 



This becomes 



f/T = cos ^/tt [1 - 2U] 



(B-11) 



with Elber’s [8] definition of the crack-closure parameter, 
U, 



U 



s — s 

max op 

S — S 

max min 



Solving Equation (B-12) for U yields 



U = I (1 - cos TTf/T) 



(B-12) 



(B-13) 



Therefore, the relation betv/een f/T and U is 



and 



1 

f/T = 2/tt sin”^ (U)^ 



(B-15) 






TABLE I. Fatigue-crack growth in 2024-T3 aluminum (from Ref. 1). 
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^Crack u liutintrd jnd |iro|^agated to a - 0.18 Inch (4.8 mm) at Sa - S ksl (34 MN/m*) to expedite testing; cycles listed are number required to propagate crack from 
0.30 Inch (7.C2 mm). 



TABLE II. Analysis of the prediction accuracy of Forman’s 
equation for small cracks in 2024-T3 aluminum. 
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TABLES III. Coefficients, di, resulting from the fit of the 
fatigue-crack growth data to the general 
second order equation' 
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Fi'gure 1. Typical fatigue-crack growth curves (from Ref. 1) 
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Figure 2. Hudson’s comparison of fatigue-crack growth rate 
with Forman's equation (From Ref. 1). 
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Mean stress = 20 ksi 
Stress amplitude = 10 ksl 
Nominal stress ratio = .3 
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Figure 3. Comparison of a fatigue-crack growth curve with the 
power-rule prediction (typical). 




Figure Comparicon of fatigue-crack growth rate with 

Forman's equation for a single load case (typical). 
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Figure 5. Constant amplitude loading with crack closure 
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Figure 6. Schematic repr^'sentation of the effect of increasingly 
negative cyclic stress ratios 
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Figure 7. Variation of the damage effective cycle fraction with cyclic 
stress ratio. 
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Figure 8. Comparison of two expressions for the crack-closure parameter 




Figure 9. Geometric Interpretation of the general second order 
polynomial fit to the fatigue-crack growth- data. 
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Figure 10. Correlation of the coefficient d-. * with the 
effective RMS gross-section stress 
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Figure 11. Correlation of the coefficient d 2 ^’ v;ith the 
effective Rr'!S gross-section sti’ess 
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Figure 12. Correlation of the coefficient with the 

effective RMS gross-section stress 
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Effective ET-S Gross-section stress, ksi 



Figure 13* Correlation of the coefficient d|j^ with 
the effective RMS gross-section stress 
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Figure Correlation of the coefficient v/lth the 

effective RMS gross-section stress 
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Figure 15. Relationship between the unstable crack propagation stress intensity 
factor and the effective RMS gross section stress 
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